Ground-state fidelity and tensor network states for quantum spin tubes 
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An efficient algorithm is developed for quantum spin tubes in the context of the tensor network representa- 
tions. It allows to efficiently compute the ground-state fidelity per lattice site, which in turn enables us to identify 
quantum critical points, at which quantum spin tubes undergo quantum phase transitions. As an illustration, we 
investigate the isosceles spin 1/2 antiferromagnetic three-leg Heisenberg tube. Our simulation results suggest 
that two Kosterlitz-Thouless transitions occur as the degree of asymmetry of the rung interaction is tuned, thus 
offering an alternative route towards a resolution to the conflicting results on this issue arising from the density 
matrix renormalization group. 
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I. INTRODUCTION 



Recently, significant progress has been made in develop- 
ing efficient numerical algorithms in the context of the ten- 
sor network (TN) representations as a consequence of our 
deeper understanding of the role of quantum entanglement 
and fidelity in characterizing critical phenomena in quantum 
lattice systems. Actually, these developments provide signif- 
icant insights into the working principle of the density ma- 
trix renormalization group (DMRG)", which in turn has led 
to some novel algorithms based on the TN representations. 
Among them are the matrix product states (MPS)' in one spa- 
tial dimension and the projected entangled-pair states (PEPS)2. 
in two or higher spatial dimensions. Especially, an infinite 
MPS (IMPS) algorithm.^ and an infinite PEPS (iPEPS) algo- 
rithm** have been developed to compute the TN representa- 
tions of ground-state wave functions for translationally invari- 
ant infinite-size quantum systems in one and two or higher 
spatial dimensions, respectively. A remarkable feature of 
these TN algorithms is that they offer an efficient way to 
compute the ground-state fidelity per lattice site'^ '^. Indeed, 
quantum fidelity allows to capture the ground-state phase dia- 
grams for various quantum lattice systems in condensed mat- 
teri^-iZ. 

There are a class of quantum lattice systems lying between 
one and two spatial dimensions, namely, quantum spin lad- 
ders and tubes. Both experiments and theories'^ have con- 
firmed that the spin 1/2 Heisenberg ladders are gapful for an 
even number of legs and gapless for an odd number of legs. 
In fact, the ground state of an odd-leg Heisenberg spin lad- 
der is a gapless spin-liquid or a Tomonaga-Luttinger liquid'''. 
Once the periodic boundary conditions are applied in the rung 
direction, spin ladders become spin tubes. Some exotic phe- 
nomena occur due to the geometric frustration for an odd-leg 
spin tube, which has been realized in some materials, e.g., for 
three-leg^" and nine-leg^' spin tubes. However, little attention 
has been paid to these systems in the context of the TN repre- 
sentations. Given the importance of both the spin ladders and 
spin tubes in the theory of strongly correlated systems, it is 
desirable to develop TN algorithms as an alternative means to 
investigate their fascinating physics. 
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FIG. 1 : (Color online) A three-leg spin tube with exchange interac- 
tion constants 7; and along the leg and rung directions, respec- 
tively. Here, an asymmetry parameter a is introduced as a tunable 
control parameter. One choice of the unit cell is highlighted in a 
dash-line box. 



This paper is a part of our effoits to develop an efficient nu- 
merical algorithm for quantum spin tubes in the context of the 
TN representations, which offers an efficient way to compute 
the ground-state fidelity per lattice site. Thus, we are able 
to locate quantum criticalities for quantum spin tubes. As 
an illustration, we investigate the isosceles spin 1/2 antifer- 
romagnetic three-leg Heisenberg tube. The simulation results 
suggest that two Kosterlitz-Thouless (KT) transitions occur as 
the degree of asymmetry of the rung interaction is tuned. 



II. TENSOR NETWORK STATES FOR QUANTUM SPIN 
TUBES 



Assume that the model Hamiltonian H describes the 
nearest-neighbor interaction and is translationally invariant 
along the leg direction. For a 3-Ieg spin tube, there are six 
sites A, B, C, D, E and F in each unit cell, as shown in Fig.[T] 
Such a choice of the unit cell is necessary to accommodate a 
ground state with spontaneously broken translational symme- 
try, as it does occur in spin tubes. Indeed, there are two differ- 
ent choices of the unit cell: one is ABCDEF, as highlighted 
in Fig.[T] the other is BADCFE. Any wave function admits a 
TN representation, which follows from attaching to each site 
a five-index tensor, as visualized in Fig. |2la) for a five-index 
tensor A^'.^^. Here, i (i = 1, 2, ■ ■ • , d) is a physical index from 
a d-dimensional local Hilbeit space. /, r, u and d denote the 
bond indices from the ;t'-dimensional auxiliary spaces, with x 
being the truncation dimension, as depicted in Fig. |2jb). In 
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FIG. 2: (Color online) A TN representation for an infinite-size three- 
leg spin tube, (a) A five-index tensor is attached to each site, with /, r, 
u and d denoting the inner indices and s the physical index, (b) A TN 
representation for a wave function, which is translationally invariant 
under two-site shifts along the leg direction, (c) An eight-index dou- 
ble tensor 05,,-,^ is formed, by contracting the physical indices, from 
the five-index tensor A^^^^ and its complex conjugate A"^,^^,^,, with 
/ = (/, /'), f = (r, r'), ft = (m, u'), and d = (d, d']. (d) The norm {ifi\t//) 
of a given wave function in the TN representation. The environment 
tensors are nothing but the left and right dominant eigenvectors of 
the transfer matrix consisting of six double tensors. Here, the trans- 
fer matrix is highlighted in the dash-line box. (e) The expectation 
value of a physical observation, such as the four-site Hamiltonian 
density, (f) The gradient d{i//\if/)/dAj*^j is equal to the norm with a 
tensor AJ'^^j absent. 



order to measure a physical observable, we define the eight- 
index double tensors, which are formed, by contracting the 
physical indices, from the five-index tensors and their com- 
plex conjugates, as shown in Fig. |2c). Fig. |2|d) shows how 
to compute the norm for a given wave function in the 
TN representation. Here, a key notion is the transfer matrix, 
which is highlighted in the dash-line box. The dominant left 
and right eigenvectors of the transfer matrix are nothing but 
the environment tensors, in the sense that it accommodates 
the effect of the remaining part of the infinite-size spin tube 
system, once the system is partitioned into two parts, with one 
of them being the transfer matrix. Notice that the norm fol- 
lows from the dominant eigenvalue Ai of the transfer matrix. 



III. UPDATING PROCEDURE FOR THE TENSOR 
NETWORK 



We need to update all the tensors in the unit cell during 
the iteration to generate the ground-state wave function for 
a given Hamiltonian. Here, we focus on how to update the 
tensor A''^^^ to explain the procedure. In fact, other tensors are 
updated in exactly the same way. 

For a given quantum state \ifr}, the energy per site, e, is writ- 
ten as 



lim 1 

L^oo L {4i\\p) 



Here, L is the total number of the lattice site. In the TN 
representation, the energy per site, e, is a functional of the 
TN tensors, which may be computed efficiently, as shown in 
Fig. He). 

Now we turn to the computation of the gradient of the en- 
ergy functional with respect to A*: 



lim ir 1 mm^) 

L™ L^(l/r|l/r) dA* 



<lA|//|tA) dim 

dA* 



], (2) 



where A* is the complex conjugate of A. Here and hereafter, 
we have omitted the tensor indices for clarity. In order to com- 
pute the energy gradient efficiently, we take advantage of the 
translational invariance of the Hamiltonian. Since there are 
two different choices of the unit cell, the contribution to the 
energy gradient involves six Hamiltonian densities acting on 
different plaquettes, i.e., Habcd, hsADC, hcoEF, hocFE, hsFAB 
and hfEBA- Then the energy gradient becomes 



with 



G,„ - 



1 d{ip\hM) <(A|/i«#> dim 



im dA* 



im^ OA* 



(3) 



(4) 



(1) 



Here, /z,„ 6 {hABCDJ''BADC^ hcDEFJ''DCFEJ''EFABJlFEBA}- A.S 

an example, we take the Hamiltonian density h^BCD to explain 
how to compute the gradient. As shown in Fig. [3] there is a 
TN representation for the gradient G^^^^, in which a hole cell 
is the cell with the tensor A* absent, and a Hamiltonian cell is 
the cell with the Hamiltonian density Iiabcd sandwiched by 
attaching to the physical indices between four five-index ten- 
sors and their complex conjugates. The expectation value of 
the Hamiltonian density /iabcd and the gradient diiff\iff} / dA* 
are shown in Fig. |2le) and (f), respectively. Notice that the 
gradient dit//\iJ/)/dA* is nothing but the norm with a tensor A* 
removed. In addition, diiff\h„,\iff}/dA*is nothing but the expec- 
tation value of h,„ with the tensor A* removed. Depending on 
the relative locations of the hole and the Hamiltonian cells, 
the TN representation for the energy gradient consists of three 
parts: (I) the hole and the Hamiltonian density posit on the 
same cell; (II) the Hamiltonian cell posits on the left hand 
side of the hole cell; and (III) the Hamiltonian cell posits on 
the right hand side of the hole cell. It should be emphasized 
that the gradient decays very fast when the Hamiltonian cell 
falls apart from the hole cell. That is, only a few Hamiltonian 
cells need to be taken into account to meet a preset precision. 
Similarly, one may compute the contribution to the gradient 
Ga from other Hamiltonian densities. 

Once one knows the energy gradient, the real and imaginary 
parts. Are and A„„, of the tensor A may be updated as follows: 
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FIG. 3: (Color online) The TN representation for the energy gradient 
G'* with respect to the tensor A*. A small circle represents a double 
tensor attached to each site, a square indicates the absence of the 
tensor A*, and a large circle is the Hamiltonian density acting on a 
plaquette in the unit cell. Here, Ai denotes the dominant eigenvalue 
of the transfer matrix. The TN for the energy gradient consists of 
three parts: (I) the hole and the Hamiltonian density posit on the 
same cell; (II) and (III) the Hamiltonian cell posits on the left and 
right hand sides of the hole cell, respectively. Here, m is the cell 
number between the hole cell and the Hamiltonian cell. Notice that 
only the contribution from the plaquette Hamiltonian density Habcd 
is shown. 



model exhibits two limiting cases: for a -0,it corresponds to 
the three-leg ladder; for o- ^ oo, it becomes a decoupled sys- 
tem consisting of a single chain and a two-leg ladder. Both of 
them are gapless Tomonaga-Luttinger liquids, with the cen- 
tral charge c - I, and possess the SU(2) symmetry. When 
a - I, the system is spontaneously dimerized, arising from 
the geometric frustration. All the spin excitations are gap- 
ful, as follows from the Lieb-Schultz-Mattis theorem^^: two 
degenerate ground states arise from the broken translational 
symmetry. Therefore, one may expect that at least one phase 
transition point occurs if the asymmetric parameter a is tuned 
from to oo. From now on, we fix 7/ and J,- to be unity, and 
focus on the effect of the asymmetry of the rung interaction. 

Previous studies for this model based on the DMRG simula- 
tions leads to somewhat controversial results: Nishimoto and 
Arikawa suggested that the spin gap vanishes as soon as the 
infinitesimally small asymmetry is introduced^^, while Sakai 
et. el}'^ argued that a finite spin gap appears in a narrow re- 
gion around the rung-symmetry line and the phase transitions 
between the gapless and gapful phases belong to the KT uni- 
versality class. 



V. NUMERICAL SIMULATION 



where G^^ and G^^ are, respectively, the real and imaginary 
parts of the five-index tensors G^, 5 is the step size for an up- 
date. During the implementation, it is tuned to be decreasing 
with the number of the iteration steps. Notice that both the 
sign and magnitude of the energy gradients G^^ and Gf^^ are 
exploited. As such, the tensor A is updated to A = + iAj,,,. 
Actually, all the six tensors in the unit cell are updated si- 
multaneously. Repeating the procedure until the ground-state 
energy per site converges, an approximate ground-state wave 
function is generated in the TN representation. 



IV. MODEL 

We consider the antiferromagnetic spin 1/2 three-leg 
Heisenberg tube, as depicted in Fig.lT] on an infinite-size lat- 
tice. The Hamiltonian takes the form 

H - //leg + //mng, (6a) 

Hieg - J' ^'J ■ (6b) 

ij 

^rung - -^r ^ ( ■ S'^Z + S;,2 ' S,^ + a S;,3 ■ S,;i ), 

(6c) 

where S,j is the spin- 1/2 Pauli operators at leg / and rung j 
(j =1,2, 3), Ji is the neighboring exchange interaction along 
the legs and Jr stands for the rung interaction. For the antifer- 
romagnetic coupling, the model is frustrated along the rung 
direction. The asymmetric parameter a controls the strength 
of the frustration. If the asymmetric parameter a is varied, the 



For a quantum system described by the Hamiltonian H(a), 
with a being the control parameter, the fidelity F{ai,a2) = 
\{^'{a2)\^{ai))\ between two ground states \^'{al)} and 
\'^(a2)} scales as F{ai,a2) ~ d(ai,a2)^. Here, d{ai,a2) 
is the ground-state fidelity per lattice site, which character- 
izes how fast changes when the thermodynamic limit is ap- 
proache d'^i'^ . It satisfies the properties inherited from the fi- 
delity F(aua2)' (i) normalization ii(a, a) = 1; (ii) symmetry 
d{ai,a2) = d{a2,aiy, and (iii) range < d(ai,a2) < 1. The 
ground-state fidelity per lattice site exhibits singularities as a 
pinch point when the control parameter a crosses a transition 
point ac- 

The TN algorithm for quantum spin tubes, as described, 
allows to efficiently compute the ground-state fidelity per lat- 
tice site. In Fig. HI the ground-state fidelity per lattice site, 
d{a[ , Q'2), is plotted for the spin 1/2 three-leg Heisenberg tube, 
with truncation dimension x - (i- Two transition points are 
identified, as two pinch points on the fidelity surface clearly 
indicates that there are two phase transitions at a^i ~ 0.95 
and ac2 ~ 1.06, respectively. Between two phase transi- 
tion points, a finite gap opens, which follows from the Lieb- 
Schultz-Mattis theorem, due to spontaneous symmetry break- 
ing of the translational symmetry, as reflected in the two-site 
reduced density matrices. Given that both gapless phases are 
Tomonaga-Luttinger liquids (see below for a detailed discus- 
sion), the phase transitions between the gapless and gapful 
phases belong to the KT universality class. Our results sup- 
port the DMRG results by Sakai et. el}'^. 

In Fig. |5] we plot the ground-state fidelity per lattice site, 
d(a\,a2), as a function of a\, for a fixed q'2, with diff'erent 
truncation dimensions Here, two reference states ^{02 - 0) 
and *I'(q'2 - 2) in the pseudo-5 f/(2)-symmetry -broken phases 
have been chosen, respectively^^. The simulation is per- 
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FIG. 4: (Color online) The ground-state fidelity per lattice site, 
d(ai,a2), for the spin 1/2 three-leg Heisenberg tube with truncation 
dimension x = 6- The fidelity surface clearly indicates that there are 
two transition points, as two pinch points occur at a^i ~ 0.95 and 
0-^2 ~ 1.06, respectively. 



allows to efficiently compute the fidelity per lattice site, which 
enables us to identify quantum criticalities for quantum spin 
tubes. Our simulation results for the isosceles spin 1/2 anti- 
ferromagnetic three-leg Heisenberg tube suggest that two KT 
transitions occur as the degree of asymmetry of the rung in- 
teraction is tuned, thus offering an alternative route towards a 
resolution to the conflicting results on this issue arising from 
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formed for a randomly chosen initial state. The algorithm au- 
tomatically produces degenerate ground states which pseudo- 
break the S U (2) symmetry in the Tomonaga-Luttinger liquid 
phases. In the pseudo-symmetry-broken phase, the ground- 
state fidelity per lattice site, d{ai,a2), exhibits different val- 
ues lying between two (extreme) degenerate ground states. In 
the symmetric phase, d{a[,a2) yields just one value. Thus, a 
bifurcation point occurs, which coincides with the transition 
point Off. Therefore, two phase transition points at a^i ~ 0.95 
and Q'c2 ~ 1.06 are identified from the bifurcation in the 
ground-state fidelity per lattice site. A smoking-gun signa- 
ture of the Tomonaga-Luttinger liquids is that such a bifurca- 
tion tends to disappear when the truncation dimension x goes 
to infinity^^, as required to keep consistent with the Mermin- 
Wagner theorem^^ that no continuous symmetry is sponta- 
neously broken in a one-dimensional quantum system. This 
implies that the degenerate ground states arise from the finite- 
ness of the truncation dimension, an artifact of the algorithm. 
However, the singularity associated with the bifurcation still 
remains (as an essential singularity as it should be for the KT 
transitions), even if the truncation dimension goes to infinity. 
Notice that the Goldstone modes survive as gapless excita- 
tions in the limiting case without truncation. 
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FIG. 5: (Color online) The ground-state fidelity per lattice site, 
d(ai,a2), as a function of ffi for a fixed a2. The truncation dimen- 
sion is fixed as = 4 and x = (>, respectively. Here, the reference 
state has been chosen: (a) 'i'(a2 = 0) and (b) *I'(q'2 = 2). 



the DMRG simulations. 
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